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We present the first measurements of tlie ^(770)", K*(892), A(1232)++, E(1385), and A(1520) 
resonances in d+Au collisions at ^s„„ = 200 GeV, reconstructed via their hadronic decay channels 
using the STAR detector at RHIC. The masses and widths of these resonances are studied as a 
function of transverse momentum (pr)- We observe that the resonance spectra follow a generalized 
scaling law with the transverse mass (mT). The (pr) of resonances in minimum bias collisions is 
compared to the (pr) of vr, K, and p. The p°/tt~, K'/K', A++/p, E(1385)/A, and A(1520)/A 
ratios in d+Au collisions are compared to the measurements in minimum bias p + p interactions, 
where we observe that both measurements are comparable. The nuclear modification factors {RdAu) 
of the p" , K* , and E* scale with the number of binary collisions (Nun) for pr > 1.2 GeV/c. 
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PACS numbers: 25.75.Dw, 25.75.-q, 13.75.Cs 



I. INTRODUCTION 

Quantum chromodynamics (QCD) predicts that 
hadronic matter at high temperatures and/or high den- 
sities undergoes a phase transition to a system of decon- 
fined partonic matter, the Quark Gluon Plasma (QGP) 
Matter under such extreme conditions can be stud- 
ied in the laboratory by colliding nuclei at very high en- 
ergies. The Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory has provided a variety 
of collision systems at different beam energies, including 
collisions of Au-|-Au, d+Au and p + p at y/s^ = 200 
GcV. 

Resonances are strongly decaying particles with life- 
times X c that are of the order of the size of the hot 
and dense medium produced in heavy-ion collisions. As 
such, the measurement of resonances in Au-l-Au collisions 
compared to p + p collisions at y/s^^ = 200 GeV has 
provided detailed information about the interaction dy- 
namics in relativistic heavy- ion collisions 0, 0, 0, S @, 0] j 
where hadronic lifetimes and interaction cross-sections 
affect the resonance yields 0, i, i, E [H E, HI • 

The in-medium effects related to the high density 
and/or high temperature of the medium can modify the 
properties of short-lived resonances, such as their masses, 
widths, and even their spectral shapes [l3,[ilj[3- Thus, 
resonances provide a unique tool for studying various 
properties of interaction dynamics in relativistic heavy- 
ion coUisions [13, [Hi- A good understanding of resonance 
production in the reference systems p + p and d-|-Au is 
useful in understanding resonance production in Au-|-Au 
collisions. Comparisons between p + p and Au-|-Au for 
the p° and K*^ mesons have been made elsewhere [1, 3] , 
where it was observed that there were modifications of 
the resonance properties (mass and width) in both sys- 
tems with respect to values in the vacuum in the absence 
of any medium effects. The measurement of masses and 
widths of resonances in d-\-Au collisions add further in- 
formation to the existing measurements. 

In addition, the regeneration of resonances and the re- 
scattering of their daughters are two competing effects 
that affect the interpretation of resonance production. 
Resonances that decay before kinetic freeze-out (vanish- 
ing elastic collisions) may not be reconstructed due to the 
re-scattering of the daughter particles. In this case, the 
resonance survival probability is relevant and depends 
on the time between chemical and kinetic freeze-outs, 
the source size, and the resonance transverse momen- 
tum (pt)- On the other hand, after chemical freeze-out 
(vanishing inelastic collisions), elastic interactions may 
increase the resonance population compensating for the 
ones that decay before kinetic freeze-out. The case of 
resonance regeneration depends on the hadronic cross- 
section of their daughters. Thus, the study of resonances 
can provide an independent probe of the time evolution 



of the source from chemical to kinetic freeze-out and de- 
tailed information on hadronic interaction at later stages. 
This has been measured in Au-|-Au and compared to p+p 
for the X*, S*, and A* 0,0] resonances. Now, with the 
addition of the rf+Au measurement we can gain insight 
into the re-scattering processes in p + p and Au+ Au col- 
lisions. 

In this paper, we present the first measurements of 
p{770)°, K*{892), A(1232)++, S;(1385), and A(1520) 
in d-l-Au collisions at y^s^ = 200 GeV at RHIC re- 
constructed via their hadronic decay channels using the 
STAR detector. The p(770)°, K*{892), A(1232)++, and 
S(1385) masses are presented as a function of px in d+Au 
collisions and the and A++ masses are compared to 
the measurements in p + p collisions. The pr spectra of 
these resonances are presented for different centralities 
in d-|-Au collisions. The (px) of resonances measured in 
minimum bias collisions are compared to the (px) of tt, 
K, and p. The p^/tt", K*/K-, A++ /p, S(1385)/A, and 
A(1520)/A ratios in d+Au and minimum bias p + p col- 
lisions arc compared. The nuclear modification factors 
RdAu of these resonances are compared to the RdAu of 
charged hadrons and the Cronin (initial state multiple 
scattering) enhancement [l^ is discussed. 



II. EXPERIMENT 

We present measurements of resonances via their 
hadronic decay channels (see Table [l| in d-l- Au collisions 

a/^jvjvL^ 200 GeV using the Time Projection Chamber 
(TPC) [201, which is the primary tracking device of the 
STAR experiment. 

A minimum bias trigger was defined by requiring that 
at least one beam-rapidity neutron impinges on the Zero 
Degree Calorimeter (ZDC) [2l| in the Au beam direc- 
tion. The measured minimum bias cross-section amounts 
to 95 ± 3% of the total c?-|-Au geometric cross-section. 
Charged particle multiplicity within the pseudorapidity 
interval —3.8 < 77 < —2.8 was measured in a forward 
TPC (FTPC) [Hi in the Au beam direction and served 
as the basis of our d-|-Au centrality tagging scheme, as 
described elsewhere [2^. The d-l-Au centrality definition 
consists of three event centrality classes; 0-20%, 20-40%, 
and 40-100% of the total d+An cross-section The 
analysis of different centralities was restricted to events 
with a primary vertex within 50 cm of the center of the 
TPC along the beam direction to ensure uniform accep- 
tance in the range studied. For the minimum bias 
events, the FTPC was not used, so events with a primary 
vertex within 100 cm were accepted, still maintaining a 
uniform acceptance in the 77 range studied. In order to 
improve the statistics in the case of the K*{8%2), events 
with a primary vertex within 75 cm were accepted for 
the centrality studies. The difference between using 50 
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TABLE I: The resonances measured in d+Au collisions, their 
corresponding hadronic decay channels, branching ratios, and 
lifetimes 



— 

Resonance 


— ;r ;:;7 ; 

Decay Channel 


■R r? 


CT 


P (770) 


^+.^- 
n n 




i.o tm 


K (892) 


i\ 7r 


fid vo/ 
~ DD. / /o 


o.y tm 


IV- + 


~ 66.7% 


4 tm 


Jf*(892)± 


KgTV 


~ 66.7% 


4 fm 


A(1232)++ 




~ 100% 


1.6 fm 


A(1232)-- 




~ 100% 


1.6 fm 


E(1385)+ 


Att+ 


~ 87% 


5.5 fm 


E(1385)- 


Att" 


~ 87% 


5.5 fm 


E(1385)- 


Att" 


- 87% 


5.0 fm 


E(1385)+ 


An+ 


~ 87% 


5.0 fm 


A(1520) 


pK- 


~ 22.5% 


12.6 fm 


A(1520) 


pK+ 


~ 22.5% 


12.6 fm 




p (GeV/c) 

FIG. 1: (Color online) dE/dx for charged particles vs. mo- 
mentum measured by the TPC in d+Au collisions. The curves 
are the Bichsel function [2^ for different particle species. 



TABLE IL The data-set for each centrality used in the anal- 
ysis of resonances in d-|-Au collisions. 



Centrality 


Number of 


Primary 


Resonance 




Events 


Vertex (cm) 




Minimum Bias 


~ 16 X 10" 


± 100 




Minimum Bias 


~ 15 X 10^ 


± 75 


K* 


Minimum Bias 


~ 14 X 10*^ 


± 70 


A* 


Minimum Bias 


~ 11.6 X 10^ 


± 50 


A++ 



and 75 cm as the primary vertex cut was taken into ac- 
count in the systematic errors. The same primary vertex 
cut was used for the minimum bias events. In order to 
improve statistics in the case of the A(1520), events with 
a primary vertex within 70 cm were accepted for cen- 
trality selected minimum bias events. A summary of the 
relevant data-sets is given in Table [III 

We also present measurements in p + p collisions of the 
A"'""'" where a minimum bias trigger was defined using co- 
incidences between two beam-beam counters that mea- 
sure the charged particle multiplicity at forward pseu- 
dorapidities (3.3 < \r]\ < 5.0). In this case, ~ 6 x 10^ 
events were used, where only events with a primary ver- 
tex within ±50 cm were accepted. 

In addition to momentum information, the TPC pro- 
vides particle identification for charged particles by mea- 
suring their ionization energy loss (dE/dx). Fig. [T]shows 
dE/dx as a function of momentum (p) measured in the 
TPC. The different bands presented in Fig. [T] represent 
Bichsel distributions folded with the experimental reso- 
lutions and correspond to different particle species [25| . 
Charged pions and kaons can be separated in momenta 
up to about 0.75 GeV/c, while (anti-)protons can be iden- 
tified for momenta of up to about 1.1 GeV/c. In Fig. 
[U the Bichsel function is used instead of the tradi- 
tional Bethe-Bloch parametrization (26j in order to im- 
prove particle identification. To quantitatively describe 
the particle identification, the variable N^y-n, which ex- 
presses energy loss in the units of the standard deviation 



of a Gaussian formed by the logarithm of truncated en- 
ergy loss, is defined (in this case for pions) as: 

Tvr 1 1 dE / u^^^mcasurcd /-, \ 

Na^ = TT T log , , . , (1) 

(^dE/dxyL-Tvc) {dE/dx},, 

where dE / dxnioa.suiod is the measured energy loss 
for a track, (dE/dx)^ is the expected mean energy 
loss for charged pions with a given momentum, and 
<^dE/dxiL'TPc) is the dE/dx resolution that depends on 
the track length in the TPC that is used in the dE/dx 
measurement. For Xtpc equal to 72 cm, corresponding 
to a 90° angle with the beam axis, the resolution is 8.1%. 
In the case of charged kaon and charged proton identifica- 
tion, similar definitions of N^-k and N^-p can be obtained. 
In order to quantitatively select on charged pions. kaons, 
and protons, specific analysis cuts, described later, are 
then applied to the variables N^„, No-k, and Na-p. 

III. PARTICLE SELECTION 

In all cases, particles and anti-particles are combined 
in order to improve statistics. In the following, the term 
stands for K*" or K*'^ , and the term K* stands 
for A'*", K* or K*^, unless otherwise specified. The 
term A"*""*" stands for A"*""*" or A , the term S* stands 
for E(1385)+, E(1385)-, S(138_5)- or E(1385)+, and the 
term A* stands for A(1520) or A(1520), unless otherwise 
specified. 

As these studied resonances decay in such short times 
that the daughters seem to originate from the interac- 
tion point, only charged pion, kaon, and proton candi- 
dates whose distance of closest approach to the primary 
interaction vertex was less than 3 cm were selected. Such 
candidate tracks are referred to as primary tracks. In or- 
der to avoid the acceptance drop in the high rj range, all 
track candidates were required to have |?7| < 0.8. For all 
candidates, in order to avoid selecting split tracks, a cut 
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FIG. 2: The reconstruction efficiency x detector acceptance 
as a function of the invariant mass for minimum bias d+Au. 
The error shown is due to the available statistics in the sim- 
ulation. 



on the ratio of the number of TPC track fit points and 
the maximum possible points was required. In addition, 
a minimum pT cut was applied to maintain reasonable 
momentum resolution. 

In the case of the a series of cuts was applied to 
the charged pion candidates in order to ensure track fit 
quality and good dE/dx resolution. A compilation of the 
cuts used in the p° analysis is given in Table Hill and the 
correction factor (reconstruction efficiency multiplied by 
the detector acceptance) as a function of invariant mass 
for a particular pT bin is depicted in Fig. [2l In general, 
the correction factor increases as a function of transverse 
momentum. The fact that the correction factor is larger 
at low values of M^^r and larger values of pr is simply 
due to kinematics. In the case of wide resonances, such as 
the p° and the A"*"*", the correction factor depends on the 
invariant mass for each pT interval that is being analyzed. 
In this case, the correction is applied as a function of the 
invariant mass for each px bin. In the case of narrow 
resonances, such as the K*, S*, and A*, the correction 
factor is dependent only on the pr bin being analyzed. 
Therefore, the correction is performed as a function of 
PT only. 

For the K* analysis, charged kaon candidates were 
selected by requiring IN^tk] < 2 while a looser cut 
\NcrT:\ < 3 was apphed to select the charged pion can- 
didates to maximize statistics for the K* analysis. Such 
cuts do not unambiguously select kaons and pions, 
but do help to reduce the background significantly. The 
background was reduced further by selecting only kaon 
candidates with p < 0.7 GeV/c. This momentum cut en- 
sures clearer identification by minimizing contamination 
from misidentified correlated pairs and thus the system- 
atic uncertainties [1]. The K*^ first undergoes a strong 
decay to produce a Kg and a charged pion hereafter la- 
belled as the K*^ daughter pion. Then, the produced Kg 
decays weakly into a tt+tt" pair with cr = 2.67 cm. The 
oppositely charged pions from the Kg decay are called 
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FIG. 3: (Color online) The K*° and if± reconstruction 
efficiency X detector acceptance as a function of pr for min- 
imum bias d-|-Au and three different centralities. 
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FIG. 4: (Color online) The A"'"^ reconstruction 
efficiency X detector acceptance as a function of the in- 
variant mass for minimum bias d-l-Au and p + p. The error 
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the K*^ granddaughter pions. The charged daughter 
pion candidates were selected from primary track sam- 
ples and the Kg candidates were selected through their 
decay topology [23, • The procedure is briefiy out- 
lined below. The granddaughter charged pion candidates 
were selected from tracks that do not originate from the 
primary collision vertex. Oppositely charged candidates 
were then paired to form neutral decay vertices. When 
the Kg candidate was paired with the daughter pion to 
reconstruct the charged K* , tracks were checked to avoid 
double counting among the three tracks used. Cuts were 
applied to the daughter and granddaughter candidates to 
ensure track fit quality and good dE/dx resolution and to 
reduce the combinatorial background in the K'g invariant 
mass distribution. All the cuts used in this K* analysis 
are summarized in Table IIIII and the K* reconstruction 
correction factors are shown in Fig. [31 

The cuts applied to the A++ and A* decay daughters 
were the same as described above for the p° and 
and their respective values are given in Table IIVI We 
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TABLE III: List of track cuts for charged kaons and charged pions and topological cuts for neutral kaons used in the and 
K* analyses in d+Au collisions. decayLength is the decay length, dcaDaughters is the distance of closest approach between 
the daughters, dcaVOPrmVx is the distance of closest approach between the reconstructed Kg momentum vector and the 
primary interaction vertex, dcaPosPrmV x is the distance of closest approach between the positively charged granddaughter 
and the primary vertex, dcaNegPrmVx is the distance of closest approach between the negatively charged granddaughter and 
the primary vertex, Mj^o is the K% invariant mass in GeV/c^, N FitPnts is the number of fit points of a track in the TPC, 
NTpcHits is the number of hits of a track in the TPC, MaxPnts is the number of maximum possible points of a track in the 
TPC, and DC A is the distance of closest approach to the primary interaction point. The Normalization Region corresponds 
to the interval in which the invariant mass and the background reference distributions are normalized. 
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as a function of pt for minimum bias d-|-Au. The error shown 
is due to the available statistics in the simulation. 
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different centralities. The error shown is due to the available 
statistics in the simulation. 



also present the A+^" measured in p + p collisions and 
the cuts and their respective values applied to the decay 
daughters are the same as the ones used in the d-|-Au 
analysis and shown in Table HVl Similarly to the K*^, 
the S*^ first undergoes a strong decay to produce a A, 
which subsequently decays weakly into n^p with a cr = 
7.89 cm. The cuts apphed to the T,*^ decay daughters 
and granddaughters are the same as mentioned for the 
K*"^ and the values are shown in Table [TVl The A++, 
A*, E*^ reconstruction correction factors arc shown in 
Figs, mini and [HI respectively. 



IV. ANALYSIS AND RESULTS 

The p° measurement was performed by calculating the 
invariant mass for each tt'^tt~ pair in an event which 
passed the cuts. The resulting invariant mass distribu- 
tion was then compared to a reference distribution calcu- 
lated from the geometric mean of the invariant mass dis- 
tributions obtained from uncorrelated like-sign pion pairs 
from the same events The 7r"''7r~ invariant mass dis- 
tribution (Af^^) and the like-sign reference distribution 
were normalized to each other between 1.5 < M^^^ < 2.5 
GeV/c2. 
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TABLE IV: List of track cuts for charged kaons, charged pions, charged protons and topological cuts for lambdas used in 
the A'^+,A*, and E** analyses in d+Au collisions, decay Length is the decay length, dcaDaughters is the distance of closest 
approach between the daughters, dcaVQPrmVx is the distance of closest approach between the reconstructed A momentum 
vector and the primary interaction vertex, dcaPoaPrmVx is the distance of closest approach between the positively charged 
granddaughter and the primary vertex, dcaNegPrmVx is the distance of closest approach between the negatively charged 
granddaughter and the primary vertex. Ma is the A invariant mass in GeV /(? , N FitPnts is the number of fit points of a track 
in the TPC, NTpcHits is the number of hits of a track in the TPC, MaxPnts is the number of maximum possible points of 
a track in the TPC, 6* is the angle in the center-of-mass of one decay particle with respect to the mother particle, and DC A 
is the distance of closest approach to the primary interaction point. The Normalization Region corresponds to the interval in 
which the invariant mass and the background reference distributions are normalized. 
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The K*, E*, and A* measurements were per- 

formed using the mixed-event technique [3], where the 
reference background distribution is built with uncor- 
rclated unlike-sign kaons and pions, protons and pions, 
lambdas and pions and protons and kaons from differ- 
ent events, respectively. The background is normalized 
over a wide kinematic range (see Tables IIIII and lIVp and 
then subtracted from the corresponding invariant mass 
distribution. 



A. Masses and Widths 

The mass and width of resonances have been of great 
interest because of their possible modification in the 
medium produced in heavy-ion collisions [l8| . It is in- 
teresting to study how the resonance masses and widths 
behave in d-|-Au collisions. 

The corresponding tt'^tt~ raw invariant mass distribu- 
tion after the like-sign background subtraction for mini- 
mum bias d+Au collisions at midrapidity (|y| < 0.5) for 
a particular px bin is shown in Fig. [T] The solid black 
line in Fig. [7| is the sum of all the well defined contri- 
butions to the Mtttt distribution (hadronic cocktail) [1]. 
The Kg was fit with a Gaussian function. The uj and 
i^*(892)° shapes were obtained from the HIJING event 
generator [29| , with the kaon being misidcntificd as a pion 



in the case of the K*°. The p°(770), the /o(980) and the 
/2(1270) were fit to a BWxPS function where BW is the 
relativistic Brcit-Wigncr function [s^l 

RW = M^^MqT 

[(Af2-Af2j2+M2r2] 

and PS is the Boltzmann factor [H, [H, [sil, ^ 



to account for phase space. Here, T is the temperature 
parameter at which the resonance is emitted [l6| and 

r = To X X (4) 

M^^ ^ [M^ - Ami) 

is the width [30| , which changes as a function of momen- 
tum due to reconstruction effects. Here, Mq and I arc 
the resonance mass and spin, respectively. The masses of 
A'g. p*^, /o; and /2 were free parameters in the fit, and the 
widths of /o and /2 were fixed according to [2^ . The un- 
corrected yields of K^, p°,LO, fo, and /2 were free param- 
eters in the fit while the K*'^ fraction was fixed according 
to the K*{892)'^ —^ttK measurement, where the amount 
of contamination was determined using a detailed simu- 
lation of the TPC response using GEANT The p°. 
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FIG. 7: (Color online) The raw 7r"*"7r~ invariant mass distri- 
bution at midrapidity after subtraction of the like-sign refer- 
ence distribution for minimum bias d-|-Au collisions with the 
hadronic cocktail fit. 



w, if*^, /o, and /2 distributions were corrected for the 
detector acceptance and efficiency determined from simu- 
lation. The signal to background ratio before subtraction 
is about 1/100. 

The cocktail fit docs not reproduce the tt^tt" raw in- 
variant mass distribution at ~600 and ~850 MeV/c^, 
respectively. This is understood to be due to other con- 
tributions to the hadronic cocktail aside from what was 
described above, e.g. the /o(600) that is not very well 
established [1^1 ■ The uj yield in the hadronic cocktail fits 
may account for some of these contributions and may 
cause the apparent decrease in the /ui ratio between 
minimum bias p + p and peripheral Au+Au interactions. 

The mass is shown as a function of pT m Fig. [H]for 
minimum bias d-t-Au interactions and 0-20%, 20-40%, 
and 40-100% of the total rf+Au cross-section. A mass 
shift of ~50 McV/c^ is observed at lowp^- The width 
was fixed at Fq = 160 MeV/c^, consistent with folding 
the p° natural width (150.9 ± 2.0 MeV/c^ ^) with the 
intrinsic resolution of the detector {Spr/PT = 0.005(1 -I- 
Pt)) [ll]. The temperature parameter used in the PS 
factor was T = 160 MeV, which was also used in the p+p 
analysis 0]. In Fig. [H only the systematic uncertainty 
for the minimum bias d-l-Au measurement is depicted 
for clarity. However, the systematic uncertainty for the 
other d+Au ccntrality measurements are similar, if not 
less than the systematic uncertainty for the minimum 
bias d-l-Au measurement. 

The p° mass at |y| < 0.5 for minimum bias and three 
different centralitics in rf-l-Au collisions at y^s^vjv — 200 
GeV increases as a function of pr and is systematically 
lower than the value reported by NA27 at CERN-LEBC- 
EHS [3^ . This experiment measured the p^ in minimum 
bias p + p collisions at ^/s = 27.5 GeV for xp > 0, where 
Xi? is the ratio between the longitudinal momentum and 
the maximum momentum of the meson. In Fig. [51 the 
shaded areas indicate the p'^ mass measured in p+p colli- 
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FIG. 8: (Color online) The mass as a function of pT at 
\y\ < 0.5 for minimum bias d-|-Au interactions and 0-20%, 
20-40%, and 40-100% of the total d-|-Au cross-section (upper 
panel). The errors shown are statistical only. The comparison 
of the p" mass as a function of pt at \y\ < 0.5 measured in 
minimum bias d-|-Au, p + p, and high multiplicity p + p 0] in- 
teractions (lower panel). The brackets indicate the systematic 
uncertainty and it is shown only for the minimum bias d-f Au 
measurement for clarity. The diamonds have been shifted to 
lower values on the abscissa by 100 MeV/c in pr for clarity. 



sions (762.6 ± 2.6 MeV/c^) by NA27 [H and the dashed 
lines represent the average of the p^ mass measured in 
e+e~ (775.6 ± 0.5 MeV/c^) [H]. The p° mass measured 
in 0-20% of the total d-|-Au cross-section is slightly lower 
than the mass measured in the most peripheral ccntral- 
ity class. The masses measured in minimum bias d-|-Au, 
P + P J E^nd high multiplicity p + p 0] interactions are 
compared in Fig. [S] The comparison shows that the p^ 
mass measured in minimum bias d-|-Au and high multi- 
plicity p + p interactions are comparable. A mass shift 
of ^^70 McV/c^ was also measured in Au+Au collisions 
3. Dynamical interactions with the surrounding mat- 
ter, interference between various tt+tt^ scattering chan- 
nels, phase space distortions due to the re-scattering of 
pions forming p*^, and Bose-Einstein correlations between 
p" decay daughters and pions in the surrounding matter 
were previously given as the possible explanations Q . It 
has been proposed [3^ that the mass shift observed in 
p + p collisions is due to tttt re-scattering, which requires 
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FIG. 9: The invariant 7v'^Tr~ mass distribution after back- 
ground subtraction for minimum bias p + p collisions mea- 
sured by NA27. For details see elsewhere [34| . The vertical 
dash-dotted line represent the average of the mass 775.8 ± 
0.5 MeV/c^ measured in e"'"e~ [2^. The vertical dashed line 
is the p° mass 762.6 ± 2.6 MeV/c^ reported by NA27 [3. 



no medium. Since one also does not expect a medium to 
be formed in d-l-Au collisions, if dynamical interactions 
are also the explanation for the mass shift, then the rc- 
scattering of the p'^ with the surrounding particles must 
exist. We also observe that the p° mass is not modified 
at high pT- 

NA27 measured the p'^ in minimum bias p + p at ^/s = 
27.5 GeV for xj? > and reported a mass of 762.6 ± 2.6 
McV/c^ 



34| . The invariant tt+tt mass distribution af- 



ter subtraction of the mixed-event reference distribution 
is shown in Fig. O The vertical dash-dotted line rep- 
resents the average of the mass 775.8 ± 0.5 MeV/c^ 
measured in e+e^ collisions [l^. The vertical dashed 
line, which accounts for the phase space, is the mass 
reported by NA27 (762.6 ± 2.6 McV/c^) As shown 
in Fig. [HI the position of the p° peak is shifted by '-^ 30 
MeV/c^ compared to the p° mass in the vacuum 775.8 ± 
0.5 MeV/c2 (H. 

NA27 obtained the p" mass by fitting the invariant 
7r+7r~ mass distribution by the (BWxPS -|- BG) func- 
tion, where in this analysis, the phase space function used 
is the same as the combinatorial background (BG). NA27 
reported a mass of 762.6 ± 2.6 MeV/c^, which is ~ 10 
MeV/c^ lower than the p° mass in the vacuum. Ideally, 
the PS factor should have accounted for the shift on the 
p" peak, and the mass obtained from the fit should have 
agreed with the p" mass in the vacuum. However, just 
like in the STAR measurement, this was not the case, 
since the phase space did not account for the mass shift 
on the position of the p" peak. 

At the CERN-LEP accelerator, OPAL, ALEPH and 
DELPHI measured the £° in inclusive e+e" reactions at 
= 90 GeV P, [13, H m. OPAL reported a shift 
in the position of the p" peak by ~ 70 MeV/c^ at low 



Xp, where Xp is the ratio between the meson and the 
beam energies, and no shift at high Xp {xp ~ 1) [H, HJ- 
OPAL also reported a shift in the position of the p^ peak 
from -10 to -30 MeV/c^, which was consistent with the 
p" measurement [HI . ALEPH reported the same shift on 
the position of p° peak as observed by OPAL [s^ . DEL- 
PHI fit the raw invariant tt+tt" mass distribution to the 
(BWxPS + BG) for Xp > 0.05 and reported a p" mass 
of 757 ± 2 McV/c^ which is 7.5 standard deviations 
below the p° mass in the vacuum (775.8 ± 0.5 MeV/c^). 
As one can see, similar to NA27, DELPHI assumed that 
the phase space was described by the background func- 
tion. Bose-Einstein correlations were used to describe the 
shift on the position of p° peak. However, high chaoticity 
parameters (A ~ 2.5) were needed [s^, [13, [l^l ■ Previous 
measurements of the p mass shift and possible explana- 
tions arc discussed elsewhere [^. The masses of the p" 
and other short-lived resonances in the vacuum are ob- 
tained only in exclusive reactions and not in inclusive 
reactions where many particles are produced. 

As previously mentioned [^, one uncertainty in the 
hadronic cocktail fit depicted in Fig. [7] is the possible 
existence of correlations of unknown origin near the p*' 
mass. An example is correlations in the invariant mass 
distribution from particles such as the /o(600) which are 
not well established [1^. The uj yield in the hadronic 
cocktail fits may account for some of these contributions. 
In order to evaluate the systematic uncertainty in the p" 
mass due to poorly known contributions in the hadronic 
cocktail, the p" mass was obtained by fitting the peak to 
the BWxPS function plus an exponential function rep- 
resenting these contributions. Using this procedure, the 
p" mass is systematically higher than the mass obtained 
from the hadronic cocktail fit. This uncertainty is the 
main contribution to the systematic uncertainties shown 
in Fig. [8|and it can be as large as ~35 McV/c^ for low 
Pt- Other contributions to the systematic errors shown 
in Fig. [5] result from uncertainty in the measurement of 
particle momenta of ~3 McV/c^ and from the hadronic 
cocktail fits themselves of ^10 MeV/c^. The systematic 
uncertainties are common to all px bins and are corre- 
lated between all ccntralities in the d+Au measurements. 

Figure [TOj depicts the mixed-event background sub- 
tracted Kn and Kgir^ invariant mass distributions for 
minimum bias (i+Au interactions at midrapidity for a 
particular pT interval of the K*'^ pT and integrated over 
the full measured px range of the K*^. The signal 
to background ratio before subtraction is 1/50 for both 
cases. The solid black line corresponds to the fit to the 
relativistic p-wave Breit-Wigner function multiplied by 
the phase space (equation [2]), with T = 160 MeV, plus a 
linear function that represents the residual background. 
This comes predominantly from correlated Kit pairs and 
correlated but misidcntified pairs. A detailed study has 
been presented previously [^. 

The K* masses and widths at \y\ < 0.5 for minimum 
bias c?+Au interactions as a function of pt are depicted 
in Fig. 111! Both mass and width were obtained by cor- 
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FIG. 10: The mixed-event background subtracted Ktt raw 
invariant mass distribution for a particular K*^ pr bin (upper 
panel) and the mixed-event background subtracted A'gvr^ raw 
invariant mass distribution integrated over the K*^ pr (lower 
panel) at \y\ < 0.5 for minimum bias d-l-Au interactions. The 
dashed lines are the linear function that describes the residual 
background. 




FIG. 11: The K* mass (upper panel) and width (lower panel) 
as a function of pr at \y\ < 0.5 for minimum bias d+Au 
collisions. In the upper panel, the solid line is the PDG if*" 
mass (0.8961 GeV/c^) M. The dashed line is the PDG K*"^ 
mass (0.8917 GeV/c^) [20]. In the lower panel, the solid line is 
the K*° and iC*± widths (0.0507 GeV/c^) [H]. The brackets 
indicate the systematic uncertainties. 



recting the K* distribution for detector acceptance and 
efficiency that was determined from a detailed simulation 
of the TPC response using GEANT [H. The mass 
increases as a function of pt and at low pt (j>t <1-1 
GeV/c) the mass is significantly smaller than previously 
reported values [1^. A similar mass shift was observed 
in minimum bias p + p collisions at ^/sNN = 200 GeV Q 
and the possible explanations are the same as described 
for the p° meson. Even though a K*^ mass shift in d+Au 
collisions has not been observed before, it is important to 
note that previous measurements were ma inly interested 
in extracting the resonance cross-section [SJ]. In addi- 
tion, we observe a mass shift at low pr of ~ 10 MeV/c^, 
while previous analyses only presented the mass inte- 
grated in Pt, xp-, or Xp. The K*^ mass is in agreement 
with previous values within errors [H]. However, this 
could be due to the limited px range covered. There is 
no significant difference between the measured K* width 
and the previous values [26j . 

The main contributions to the systematic uncertainties 
on the K* mass and width were evaluated as a function 



of Pt using a different residual background function (sec- 
ond order polynomial), different fitting functions to the 
K* invariant mass (non-relativistic BW, relativistic BW 
without phase-space factor), and different slope param- 
eters in the BWxPS function (140 MeV and 180 MeV). 
In addition, the mass and the width were obtained sep- 
arately for K*°, 'K*°, K*^, and K*~ . The systematic 
uncertainty due to detector effects was also accounted 
for. The systematic uncertainty can be as large as ^6.5 
(9.5) MeV/c^ and -25 (30) MeV/c^ for the K*°{K*^) 
mass and width, respectively. 

The pTT raw invariant mass distributions after the 
mixed-event background subtraction for minimum bias 
rf-|-Au and p + p interactions at midrapidity for a par- 
ticular Pt bin are shown in Fig. [121 Before background 
subtractions, the signal to background ratios are 1/50 
and 1/30 for minimum bias c?-|-Au and p + p interactions, 
respectively. The solid black line corresponds to the fit 
to a relativistic p-wave Breit-Wigner function multiplied 
by the phase space, with T = 160 MeV, plus a Gaussian 
function that represents the residual background indi- 
cated by a dashed line. In this case, the normalization 
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FIG. 12: The mixed-event background subtracted pn raw in- 
variant mass distribution for a particular pr bin at |y| < 0.5 
for minimum bias d-l-Au collisions (upper panel) and for min- 
imum bias p + p collisions (lower panel). The dashed lines are 
the linear function that describes the residual background. 



A"^"^ Minimum Bias d+Au 




0.2 0.4 0.6 0.8 



1.2 1.4 1.6 

(GeV/c) 



FIG. 13: The A"*"^ mass (upper panel) and width (lower 
panel) as a function of at \y\ < 0.5 for minimum bias 
d-l-Au collisions. The solid lines are the PDG A"*"^ mass 
(1.232 GeV/c^) and width (0.120 GeV/c^) The brackets 
indicate the systematic uncertainties. 



factor used to subtract the combinatorial background was 
changed until the best x^/ndf was achieved. Similar to 
the analysis Q , the uncorrected yield of the A++ was 
a free parameter in the fit and the distribution was 
corrected for the detector acceptance and efficiency de- 
termined from a detailed simulation of the TPC response 
using GEANT [s^ . The relativistic p-wave Breit- Wigner 
function multiplied by the phase space is the same as 
equation [21 However, in the case of the A+"*", the mass 
dependent width is given by: 



FoA/o fc(AV)3F(A,,fc(A'V))2 



JJTT 



k{MofF{K,k{Ah)f 



(5) 



where F(A^, kcu) is the form factor used to fit the tt — N 
scattering phase-shift with A^ = 290 MeV/c^ [i^l, and 



ml) 



2 ^2 



(6) 



The A++ mass and width at |y| < 0.5 for minimum 
bias d-l-Au interactions as a function of are depicted 
in Fig. 1131 The A^+ mass is significantly smaller than 
the values previously reported, though the width is in 



agreement within errors |2a | . Possible explanations for a 
A"*""*" mass shift are the same as for the p° 0]. 

The A++ mass and width at midrapidity for minimum 
bias p + p collisions as a function of pt are shown in Fig. 
1141 The analysis procedure in minimum bias p -I- p is 
the same as in d-|-Au collisions. Similarly to the d-|-Au 
measurement, the A"''^ mass is significantly smaller than 
the values in [26[ and the same possible explanations ap- 
ply. The A"*""*" width is in agreement with previous values 
within errors [1^ . 

In the case of the A^"*" mass and width, the main con- 
tributions to the systematic uncertainties were calculated 
as a function of pT by using different residual background 
functions (first and second-order polynomial) and differ- 
ent slope parameters in the BWxPS function. The mass 
and the width were also obtained separately for A++ and 

A . The contribution due to the uncertainty in the 
measurement of particle momenta is ~5 MeV/c^. The 
systematic uncertainty can be as large as ^20 MeV/c^ 
and ~30 MeV/c^ for the A"*"^ mass and width, respec- 
tively. In p + p collisions, the systematic uncertainty on 
the mass and the width was evaluated similarly to the 
measurement in d-|-Au collisions. The systematic uncer- 
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FIG. 14: The A"'"^ mass (upper panel) and width (lower 
panel) as a function of pr at midrapidity for minimum bias 
p + p interactions. The solid lines are the PDG A^""" mass 
(1.232 GeV/c^) and width (0.120 GeV/c^) The brackets 
indicate the systematic uncertainties. 



tainty is ~10 MeV/c^ and --15 MeV/c^ for the A++ mass 
and width, respectively. 

The mid-rapidity and pr integrated Att (S*) and pK 
(A*) raw invariant mass distributions, after the mixed- 
event background subtraction, from minimum bias d-l-Au 
collisions are depicted in Fig. [15] and Fig. [1^1 respec- 
tively. The signal to background ratio is 1/14 for the 
E* and 1/24 for the A* before mixed-event background 

subtraction. Since the S~ and the have the same 
final state as the E* and E , the Att invariant mass 
distribution is fitted to a Gaussian combined with a non- 
relativistic Breit-Wigner function (SBW) 



SEW 



To 



(Ma. - AIo)^ + r2/4 



(7) 



In the case of the A*, the signal is fitted to a non- 
relativistic Breit-Wigner function combined with a linear 
function that describes the residual background 01- 

The fit parameters corresponding to the E* mass and 
width in the integrated pr interval (0.25 < pr < 3.5 
GeV/c) are 1.376 ±0.002 ±0.007 GeV/c^ and 48±2±10 
MeV/c^, respectively. Both the measured width and the 
mass, within their uncertainty, are in agreement with the 
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FIG. 15: The mixed-event background subtracted Att raw 
invariant mass distribution integrated over the E* pr at \y\ < 
0.75 for minimum bias d-l-Au collisions. 
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FIG. 16: The mixed-event background subtracted pK raw 
invariant mass distribution integrated over the A* pT at \y\ < 
0.5 for minimum bias d-l-Au collisions. The dashed line is the 
linear function that describes the residual background. 
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FIG. 17: The S* mass as a function of pt at \y\ < 0.75 for 
minimum bias d-|-Au collisions. The solid line is the PDG 
E* mass average between the masses of E*"^ and E*^ (1.3850 
GeV/c^) [2^. The brackets indicate the systematic uncer- 
tainties. 
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lution of the detector (33[ of 6 MeV and the momentum 
dependent mass shifts in the data, which are in the sta- 
tistical and systematical limits. The measured width in 
each momentum bin is consistent with folding the A* nat- 
ural width of 15.6 ± 1.0 MeV/c^ [2^ with the detector 
resolution. The systematic uncertainties arc due to the 
residual background, the range used for the normaliza- 
tion and for the fit to the signal, and different bin widths. 

In d-f-Au collisions, we observe modifications of the 
mass and decay width of short-lived resonances that 
might be due to dynamical interactions with the sur- 
rounding matter, interference, phase space, and Bose- 
Einstein correlations @. 



FIG. 18: The A* mass as a function of pr at \y\ < 0.5 for 
minimum bias d-|-Au collisions. The solid line is the PDG 
A* mass (1.5195 GeV/c^) [1^. The brackets indicate the 
systematic uncertainties. 
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FIG. 19: The A* width as a function of pr at |y| < 0.5 for 
minimum bias d-|-Au collisions. The solid line is the PDG 
A* width (0.0156 GeV/c^) [H]. The brackets indicate the 
systematic uncertainties. Our simulations, including both the 
detector resolution and kinetic cuts, show a width of 0.022 
GeV/c^ 



PDG values of the E* [26[ . The systematic errors include 
the uncertainty due to choice of bin size, the normaliza- 
tion of the mixed event background, the variations in the 
fit range and the selections of event and tracks. It is pos- 
sible to further study the pr dependence of the E* mass 
when the width is fixed to the PDG value (37.6 ± 1.1 
MeV/c^) [1^ and the mass is a free parameter in the 
Breit-Wigner function. Figure [T7] shows the px depen- 
dence for the E* mass from the fit function. There is a 
small difference in the mass for low pT E* compared to 
the PDG value. 

The results for the A* mass and width are shown in 
Figs. [TS] and [HI respectively. The A* mass obtained 
from the data is 1515.0 ± 1.2± 3 MeV/c^, consistent 
with the A* natural mass of 1519.5 ± 1.0 MeV/c^ 
within errors. The width of the pT integrated spectrum 
is 40 ± 5 ± 10 MeV/c^ which includes the intrinsic reso- 



B. Spectra 

In p+p collisions at RHIC, a shape difference in the pT 
spectra of mesons and baryons for non-resonant particles 
in the interval 2 < pr < 6 GeV/c at y/s^ = 200 GeV 
was observed [4l|. In order to verify if such effect is 
observed in d-|-Au collisions for resonances, their spectra 
are studied. 

The uncorrected yields obtained in each pT bins were 
corrected for the detector acceptance and efficiency de- 
termined from a detailed simulation of the TPC response 
using GEANT [3^ . The yields were also corrected for the 
corresponding branching ratios listed in Tabic HI to ac- 
count for the fact that wc only measure certain decay 
modes. 

The p° and the {K* + K*)/2 corrected invariant yields 
[(fN/{2TTpTdpTdy)] at \y\ < 0.5 as a function of pr for 
minimum bias d-f-Au interactions are shown in Fig. [5D] 
and Fig. [2T1 respectively. A Levy function [3] 



1 d^N 



dN 



(n-l)(n-2) 



2TrpT dydpT dy 2'KnT{nT - 



(1 



\rp¥ 



moin - 2)) 
I- — Too s 



nT 



(8) 



was used to extract the p° and K* yields per unit of 
rapidity around midrapidity. In the limit of low pr, the 
Levy function is an exponential function and a power law 
in the limit of high px- 

The (A+++A"")/2, (E*--hE*"^)/2, and {A*+T)/2 
corrected invariant yields at |?/| < 0.5 as a function of pT 
arc shown in Fig. [22l Fig. [23l and Fig. [24l respectively. 
Figure [22] also depicts the (A++-I-A )/2 corrected in- 
variant yield for minimum bias p+p. Since the pt region 
is limited to low pT, we use an exponential function [j| 



1 



d^N 



dN 



1 



2nmT dydniT dy 



exp( 



27rr(mo 

-(tot - 



-T) 
Too) 



T 



(9) 



to extract the A++, E*, and A* yields per unit of rapidity 
around midrapidity. Due to limited statistics, only the 
A* yield in minimum bias d-|-Au collisions was measured. 
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TABLE V; The p° and {K' +K')/2 dN/dy, T, and n a.t \y\ < 0.5 measured in minimum bias d+Au collisions and 0-20%, 
20-40%, and 40-100% of the total d+Au cross-section. The first error is statistical; the second is systematic. 

Resonance Centrality dN/dy T (MeV) n 

]P Minimum Bias d+Au 0.812 ± 0.004 ± 0.085 231.7 ± 1.6 ± 35 11.1 ± 0.2 

0-20% 1.169 ± 0.014 ± 0.17 245 ± 6 ± 52 13.6 ± 1.7 

20-40% 0.958 ± 0.011 ± 0.14 230 ± 4 ± 44 11.4 ± 0.6 

40-100% 0.607 ± 0.005 ± 0.13 212 ± 3 ± 36 10.7 ± 0.4 

K* Minimum Bias d+Au 0.161 ± 0.002 ± 0.027 286 ± 7 ± 44 10.4 ± 0.1 

0-20% 0.294 ± 0.009 ± 0.051 316 ± 22 ± 53 12.8 ± 0.4 

20-40% 0.204 ± 0.005 ± 0.037 306 ± 17 ± 50 12.5 ± 0.3 

40-100% 0.108 ± 0.002 ± 0.018 232 ± 7 ± 39 7.3 ± 0.6 



TABLE VL The (A++ -I- A )/2, (E* -I- E*)/2, and (A* + A*)/2 dN/dy and T at \y\ < 0.5 measured in minimum bias d-|-Au 
collisions and 0-20%, 20-40%, and 40-100% of the total d+Au cross-section. In the case of the A"*"^, the results from the 

measurements in minimum bias p + p are also shown. The first error is statistical; the second is systematic. 

Resonance Centrality dN/dy T (MeV) 

A++ Minimum Bias d+Au 0.0823 ± 0.0012 ± 0.0099 284 ± 7 ± 45 

0-20% 0.177 ± 0.005 ± 0.021 328 ±17 + 52 

20-40% 0.116 ± 0.003 + 0.014 303 + 14+ 48 

40-100% 0.0529 ± 0.0008 + 0.0063 290 + 9 + 46 

Minimum Bias p + p 0.0139 ± 0.0008 ± 0.0050 216 ±13 + 86 

E* Minimum Bias d+Au 0.0319 ± 0.0011 ± 0.0041 387 ± 11 ± 28 

0-20% 0.068 ± 0.006 ± 0.011 473 ± 39 ± 40 

20-40% 0.040 ± 0.004 ± 0.007 420 ± 36 ± 40 

40-100% 0.018 ± 0.002 ± 0.005 428 ± 36 ± 40 

A* Minimum Bias d±Au 0.0149 ± 0.0014 ± 0.0022 392 ± 75 ± 39 




and 40-100% of the total d+Au cross-section. The lines are invariant yields as a function of pr at \y\ < 0.5 for minimum 
fits to a Levy function (equation H). The errors are statistical ^ias d+Au collisions and three different centralities. The lines 
only and smaller than the symbols. ^'"^ ^° ^ Levy function (equation H) . The errors are sta- 

tistical only and smaller than the symbols. 



The extracted dN/dy, T, and n for the p° and the 
K* are listed in Table El In the case of the A++, S*, 
and A*, the corresponding dN/dy and T are listed in 
Table IVll One contribution to the systematic uncertain- 
ties quoted in Tables |V] and IVII is due to the tracking 
efficiency (~8%), which is common to all measurements. 

In the case of the , the normalization between 



the A/tttt and the like-sign reference distributions is the 
largest contribution to the systematic uncertainty to the 
yield and the inverse slope (T) and it can be as large as 
^20%. If the invariant yield is obtained for the case 
that the p° width is a free parameter in the hadronic 
cocktail, the invariant yield increases by 22% from those 
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FIG. 22: The ( A++ + A )/2 invariant yields as a function of 
Pt at \y\ < 0.5 for minimum bias p+p, d+Au interactions and 
0-20%, 20-40%, and 40-100% of the total d+Au cross-section. 
The lines are fits to an exponential function (equation|9}. The 
errors are statistical only and smaller than the symbols for the 
spectra measured in d+Au. In the p + p measurement, the 
errors shown also include the systematic uncertainties. 
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FIG. 23: The E* invariant yields as a function of pt at \y\ < 
0.75 for minimum bias d+Au collisions and three different cen- 
tralities. The lines are fits to an exponential function (equa- 
tion [Q)). 



shown in Tabic [V] In the other extreme, if the invariant 
yield is obtained by assuming an exponential background, 
the yields decrease by 45%. 

The contributions to the systematic uncertainty on 
the K* and A++ yields and T measured in d+Au colli- 
sions were obtained by comparing different BW functions 
(relativistic and non-relativistic), using different resid- 
ual background functions (first and second-order polyno- 
mial), different functions to fit the spectra (exponential 
and power-law), and different slope parameters in the 
BWxPS function (140 MeV and 180 MeV). In addition, 

the yields and T were obtained separately for K , 
K*+, K*-, A++, and A . The effect of opening the 
primary vertex from 50 cm to 75 cm in the case of the 
yields obtained for different centralitics in c?+Au coUi- 
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FIG. 24: The A* invariant yields as a function of px at \y\ < 
0.5 for minimum bias d+Au collisions. The line is a fit to an 
exponential function (equation 
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FIG. 25: (Color online) Proton and pion spectra [5^ plotted 
together with the rescaled minimum bias d+Au spectra of 
K* , A+^, and E*, where the transverse mass scaling is 
observed. The inset plot is a zoom-in of the region 1 < mr < 
2 GeV/c^ on a linear scale. 



sions was also taken into account. The systematic uncer- 
tainty on both yields and T is ~ 20% for the K* . In the 
case of the A"'"'", the systematic uncertainties are 12% 
and 17%, respectively. 

In minimum bias p + p collisions, the main contribu- 
tions to the A^-*" yield and T systematic uncertainty was 
estimated from the invariant yields as a function of pt 
by increasing the normalization between the Mp^ and 
the mixed-event reference distributions until the fit to 
the A^'--*- signal is not reasonable. This procedure is then 
repeated by decreasing the normalization. During this 
procedure, the width was fixed to 110 MeV/c^ [i^ . 

The number of partons (primarily gluons) in a nucleus 
grows very rapidly at very high energies. If the occu- 
pation number of these partons is large, they may satu- 
rate and form a novel state of matter called a color glass 
condensate (CGC). This CGC has a bulk scale which is 
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FIG. 26: (Color online) The p° and if* (pt) as a function 
of (dNch/drj) compared to that of n~ , K~ , and p for mini- 
mum bias p + p, minimum bias d+Au, and 0-20%, 20-40%, 
and 40-100% of the total d-f An cross-section [sH]- The errors 
shown are the quadratic sum of the statistical and systematic 
uncertainties. 
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FIG. 27: (Color online) The A++, E*, and A* (pr) as a 
function of (dNch/drj) compared to that of , , and p 
for minimum bias p + p, minimum bias d+Au, and 0-20%, 
20-40%, and 40-100% of the total d+Au cross-section [H]. 
The errors shown are the quadratic sum of the statistical and 
systematic uncertainties. 



the typical intrinsic transverse momentum of the satu- 
rated gluons in the nucleus. The CGC can be probed in 
deep inelastic scattering (isl . [44 1, i n photo-production in 
peripheral heavy-ion collisions 451 . in p(d)-fA collisions 
[43,|43| and in heavy-ion collisions [4^, |49|, [50| . Figure [25] 
shows that the transverse mass (mr) spectra of identified 
hadrons follow a generalized scaling law in c?-|-Au colli- 
sions between 1 < tot < 2 GeV/c^. Even though this 
scaling behavior is motivated by the idea of a saturation 
of the gluon density, the identified particle spectra mea- 
sured in p -I- p collisions at ISR [5l|, [H, [11| , at SppS [sj , 
and at RHIC [4l| energies have also been shown to follow 
a generalized scaling law in transverse mass. More theo- 
retical work is needed in order to explain the similarities 
between p + p and d+An collisions. 

It is interesting to notice that for resonances in d+Au 
collisions in the pt region measured, we do not observe 
the shape difference of the pT spectra observed for mesons 
and baryons in p + p collisions at RHIC [4l| for non- 
resonant particles in the interval 2 < pT < 6 GeV/c 
at the same beam energy. This baryon-meson effect ob- 
served inp+p collisions was argued to be a simple reflec- 
tion of the underlying dynamics of the collision in that 
meson production from fragmentation requires only a 
(quark, anti-quark) pair while baryon production requires 
a (di-quark,anti-di-quark) pair. 



The p". A'*, A++, A*, and S* (pr) as a function of 
dNch/drj compared to that of tt~ , K~, and p for min- 
imum bias o?-|-Au [s^ are depicted in Figs. [2SI and [571 
While the (px) of these hadrons are independent of cen- 
trality, as expected, the (pr) is strongly dependent on 
the mass of the particle. 

We can compare the spectra shape among particles 
for different systems by comparing their {pr}- Figure [28l 
shows the (px) of various particles for different systems, 
minimum bias p + p, d + Au, and central Au+Au colli- 
sions. Even though there is a strong mass dependence, 
the (pt) of these particles do not appear to strongly de- 
pend on the collision system, with the exception of the p. 
However, the (pt) of particles measured in rf+Au colli- 
sions lie between the (px) measured in p + p and Au+Au 
collisions, indicating a hardening of the spectra from p+p 
through d-|-Au to Au-|-Au collisions. 

The main contributions to the systematic uncertainties 
quoted in Table [VTI] are due to tracking efficiency (~8%) 
and different fit functions used to fit the pt spectra. In 
the case of the p*^, in addition there was the normalization 
between the tt~^tt~ invariant mass distribution and the 
like-sign reference distributions (~5%). 



D. Particle Ratios 



C. Mean Transverse Momenta (pr) 

The averaged transverse momentum {{pt)) provides 
information on the shape of the particle spectra. At a 
given mass, the larger the (pt), the harder the spectra 
are. The resonance (pt) were calculated from the fit 
parameters depicted in Table [Vj and Table IVII and are 
listed in Table [Vni 



It has been previously shown that the ratio of yields 
of resonances to the yields of stable particles can effec- 
tively probe the dynamics of relativistic heavy-ion colli- 
sions [2, [1, [3j [3|- The ratios of yields of resonances to 
stable particles with similar quark content but different 
spin and masses are given in Tabl e IVIIII The values of 
TT, K, and p were taken from [sl, l57j . ^Figures [211 [SOI 
[5T1 [5^ and [331 show the ratios of resonances to their cor- 
responding stable particles as a function of the charged 
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TABLE VII: The K* , A++, E*, and A* (pr) in minimum bias d+Au collisions and 0-20%, 20-40%, and 40-100% of the 
total d+Au cross-section. The {pt) measured in minimum bias p + p is also listed. The first error is statistical; the second is 
systematic. 



Resonance 


Centrality 


(pt) (GeV/c) 




Minimum Bias d-|-Au 


0.808 ± 0.050 ± 0.086 




0-20% 


0.815 ± 0.020 ± 0.083 




20-40% 


0.805 ± 0.015 ± 0.082 




40-100% 


0.764 ± 0.009 ± 0.081 




Minimum Bias p + p 


0.616 ± 0.002 ± 0.062 


K* 


Minimum Bias d+Au 


0.96 ± 0.02 ± 0.16 




0-20% 


1.00 ± 0.07 ± 0.17 




20-40% 


0.98 ± 0.05 ± 0.17 




40-100% 


0.90 ± 0.03 ± 0.15 




Minimum Bias p -I- p 


0.81 ± 0.02± 0.14 


A++ 


Minimum Bias d-|-Au 


0.89 ± 0.02 ± 0.14 




0-20% 


0.98 ± 0.05 ± 0.16 




20-40% 


0.92 ± 0.04 ± 0.18 




4U-iUU/o 


n (in -1- n oq i ri i a 
U.9U ± U.Uo ± U.i4 




Minimum Bias p -I- p 


0.63 ± 0.04 ± 0.13 


E* 


Minimum Bias d-f Au 


1.13 ± 0.03 ± 0.08 




0-20% 


1.33 ± 0.06 ± 0.10 




20-40% 


1.20 ± 0.07 ± 0.10 




40-100% 


1.22 ± 0.07 ± 0.10 




Minimum Bias p + p 


1.015 ± 0.015 ± 0.07 


A* 


Minimum Bias d4-Au 


1.17 ± 0.15 ± 0.12 




Minimum Bias p + p 


1.08 ± 0.09 ± 0.05 



1.6 

> 

a 1.4 



cT 1.2 

V 

1 

0.8 
0.6 
0.4 



0.2 




A Central Au+Au 
• MinBias d+Au 
▼ MinBias p+p 



K" 




p° K*J% AA^ 

P , 



Z A Q" 



0.2 0.4 0.6 0.8 1 



1.2 1.4 1.6 1.? 

Mass (GeV/c2) 



FIG. 28: (Color online) The (pt) of various particles for differ- 
ent systems, minimum bias p-f p, d-f Au, and central Au-f Au 
collisions. The errors shown are the quadratic sum of the 
statistical and systematic uncertainties. 



particle multiplicity {Nch) in p + d-f Au and Au-f Au 
collisions. We observe that the p"/7r^, A++/p, and E*/A 
ratios are independent of multiplicity while the K* / K~ 
and A* /A ratios seem to decrease. The resonance abun- 
dance could be affected by mass shifts due to phase space 
(exp(— to/T)) in different collision systems. 

The resonance ratios normalized by their value mea- 
sured va. p + p collisions at the same ^/s are plotted in 
Fig. [M] The decrease of the resonance ratios of K* /K~ 
and A* /A from p + p to Au-fAu collisions has been ex- 



plained by an extended lifetime of the hadronic phase 
where the re-scattering of the decay particles dominates 
over resonance regeneration 0, 0, 0, B B [111 ■ As the 
K* /K~ and the A*/A ratios are similar for p -\- p and 
d-f Au collisions, this would suggest the absence of an ex- 
tended hadronic medium in d-f Au collisions. The p^/tt", 
and E*/A ratios in d-f Au collisions are in agree- 
ment with their ratios measured in p-f p collisions. These 
resonance ratios do not show any suppression from p + p 
to Au-fAu collisions either, hence they are not sensitive 
to the lifetime of the hadronic medium, presumably due 
to their large regeneration cross-section. 

The p*^/7r~ ratio is independent of centrality up to the 
40-80% of the inelastic hadronic Au+Au cross-section 
and it is of the same order as the corresponding p + p 
measurement. In p + p collisions, it has been proposed 
that the mass shift is due to tttt re-scattering, even in 
the absence of a medium [sHi- If this is the case, tt+tt^ 
re-scattering might regenerate the In addition, one 
of the decay daughters might also re-scatter with other 
hadrons preventing the to be measured. Therefore, 
these two processes compete with (and balance) each 
other. 



E. Nuclear Modification Factor 



The nuclear modification factor {RdAu) is defined as 



RdAuiPr) 



d?NdAuldydpj 



(10) 
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TABLE VIII: The (f j-R^ , K* / , A++/p, E7A, and A7A ratios in minimum bias p + p [1,11,01, d+Au, and 0-20%, 20-40%, 
and 40-100% of the total d+Au cross-section. The first error is statistical; the second is systematic. 

Centrality p° K* / A++/p S'/A AyA 

Min. Bias d+Au 0.175 ± 0.004 ± 0.054 0.28 ± 0.01 ± 0.03 0.185 ± 0.005 ± 0.028 0.23 ± 0.03 0.106 ± 0.024 
0-20% 0.139 ± 0.014 ± 0.036 0.29 ± 0.01 ± 0.03 0.206 ± 0.006 ± 0.028 0.24 ± 0.04 

20-40% 0.158 ± 0.011 ± 0.056 0.29 ± 0.01 ± 0.03 0.192 ± 0.005 ± 0.028 0.21 ± 0.04 

40-100% 0.211 ± 0.005 ± 0.068 0.34 ± 0.01 ± 0.04 0.203 ± 0.006 ± 0.028 0.23 ± 0.06 

Min. Bias p + p 0.183 ± 0.001 ± 0.027 0.35 ± 0.01 ± 0.05 0.132 ± 0.002 ± 0.049 0.029 ± 0.047 0.092 ± 0.026 
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FIG. 29: (Color online) The p^/tt^ ratios in p + p, various 
centralities in d-|-Au, and in peripheral Au -\- Au collisions as 
a function of dNch/drj. The errors shown are the quadratic 
sum of the statistical and systematic uncertainties. 
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FIG. 31: (Color online) The A++/p ratios in p + p and vari- 
ous centralities in d+Au collisions as a function of dNch/drj. 
The errors shown are the quadratic sum of the statistical and 
systematic uncertainties. 
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FIG. 30: (Color online) The K* /K~ ratios in p+p and various 
centralities in d+Au and Au+Au collisions as a function of 
dNch/dr]. The errors shown are the quadratic sum of the 
statistical and systematic uncertainties. 



where cPNdAu/dydpT is the yield of the produced parti- 
cles in minimum bias d+Au collisions, (Pupp/ dydpT is the 
inclusive cross-section in p+p collisions, {Nun) is the av- 
erage number of binary nuclcon-nucleon (NN) collisions 
per event, and {Ni,i„) / a^p'^'' is the nuclear overlap TA{b) 
[ll,[5i,[5^. The value of ct™^' is 42 mb. 

The enhancement observed in RdAu for high pT and 
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FIG. 32: (Color online) The E*/A ratios in p + p, various 
centralities in d+Au, and in central Au+Au collisions as a 
function of dNch/drj. The errors shown are the quadratic 
sum of the statistical and systematic uncertainties. 



mid-rapidity, known as the Cronin effect [l9|, is gener- 
ally attributed to the influence of multiple parton scat- 
tering through matter prior to the hard scattering that 
produces the observed high-p^ hadron (60j . Therefore, 
the nuclear modification factor (RdAu) can be used to 
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FIG. 33: (Color online) The A*/A ratios in p + p, minimum 
bias d+Au, and two different centralities in Au+Au collisions 
as a function of dNch/dr^- The errors shown are the quadratic 
sum of the statistical and systematic uncertainties. 
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FIG. 35: The RdAu compared to the charged hadrons 
RdAu- The shaded box is the error on the overall normal- 
ization and the black box is the error on Ntin- 



< 
a: 



2.2 
2 
1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



(h +\\)I2 

n'+n* 

• K*0 

O K*± 



O 



■■xr 



0.5 



1.5 



2.5 



3 5 4 
p.|.(GeV/c) 



FIG. 34: (Color online) The resonance ratios normalized by 
their ratio measured in p + p collisions at the same beam 
energy as a function of dNch/dr). The errors shown are the 
quadratic sum of the statistical and systematic uncertainties. 



FIG. 36: The K* RdAu compared to the charged hadrons 
RdAu- The shaded box is the error on the overall normaliza- 
tion and the black box is the error on Nbin- 



study the effects of matter on particle production. 

The RdAu for p°, K* , and E* are shown in Figs. [351 
[36l and [37] together with the RdAu of charged hadrons 
and charged pions. The K* and S* RdAu are lower than 
unity at low pT and consistent with the RdAu of charged 
hadrons and charged pions. The RdAu of the K* , 
and S* scale with Nbm for pr > 1-2 GeV/c taking into 
account the uncertainties in the normalization. We also 
observe that the RdAu for px > 1-5 GeV/c is sup- 
pressed compared to the charged hadrons and charged 
pions RdAu- The A"*""*" RdAu is not shown due to the 
small pt range covered and the large uncertainties in the 
measurement. 

More information may be obtained from the RdAu mea- 
surement if it is extended to higher pT- In STAR this 
will be possible with the installation of the barrel Time- 
of-Fhght (TOP) detector. The TOF will provide essen- 
tial particle identification by, for instance, increasing the 
percentage of kaon and protons for which particle iden- 



tification is possible to more than 95% of all those pro- 
duced within the acceptance of the TOF barrel (jTyj < 
1.0). The improvement in particle identification will al- 
low a decrease in the signal to background ratios for the 
resonance measurements. 



V. SUMMARY 

Measurements of p(770)", /\*(892), A(1232)++, 
E(1385), and A(1520) in 200 GeV d+Au col- 

lisions reconstructed via their hadronic decay channels 
using the STAR detector are presented. 

The masses of the pO, K*, A++, E(1385), and A(1520) 
are measured for minimum bias and three different cen- 
tralities in (i+Au collisions at ^/snn = 200 GeV. We 
observe a p° mass shift of ~50 MeV/c^ at low pr- In 
addition, the /9° mass measured in 0-20% of the total 
d+Au cross-section is slightly lower than the mass mea- 
sured in the most peripheral centrality class. We also ob- 
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FIG. 37: The E* RdAu compared to the charged hadrons 
RdAu- The shaded box is the error on the overall normaliza- 
tion and the black box is the error on Ntin- 



serve that the mass measured in minimum bias d+Au 
and high multiplicity p + p interactions are comparable. 
The K*''^ and S* masses at low px [pt <1-1 GeV/c) are 
smaller than previously measured values [1^ by up to 
~10 McV/c^. A similar mass shift for the K*^ is ob- 
served in minimum bias p + p collisions at ^Jsmn = 200 
GeV The K*"^ mass and the K* width are in a gree - 
ment with previously reported values within errors [26| . 
The A++ mass is shifted by ^40 MeV/c^ in minimum 
bias p + p and ^--^50 MeV/c^ in minimum bias d+Au col- 
lisions. In contrast to the no pr dependence is ob- 
served. Similar mass and/or width modifications with re- 
spect to those observed in e-|-e- collisions are observed for 
these resonances in p+p and Au-|-Au collisions. The pos- 
sible explanations for the apparent modification of the 
meson properties are interference between various Tr~^T:^ 
scattering channels, phase space distortions due to the 
re-scattering of pions forming p'^ , and Bose- Einstein cor- 
relations between decay daughters and pions in the 
surrounding matter Q. All these explanations require 
an interaction, which implies a medium such as the one 
formed in A-l-A collisions. However, the p° mass shift 
measured in p+p collisions 0] can be described without a 
medium [s^. The A+"'" width and the A* mass and width 
measured are in agreement with previous measurements 

The transverse mass spectra follows a generalized scal- 
ing law between 1 and 2 GeV/c^. However, in Au-t-Au 
collisions at RHIC at ^snn = 200 and 62.4 GeV a gen- 
eralized scaling law is not observed [sH , possibly due to 
additional physics effects such as flow, coalescence and 
energy loss. Even though the scaling behavior in d+Au 
collisions is motivated by the idea of a saturation of the 
giuon density, the identified particle spectra measured in 
p+p have also been shown to follow a generalized scaling 
law in transverse mass. More theoretical work is needed 
in order to explain the similarities between p + p and 
d-|-Au collisions. The resonances in d+Au collisions in 
the measured px region do not show the shape difference 



of the Pt spectra observed for mesons and baryons in 
p + p collisions at RHIC for non-resonant particles in the 
interval 2 < pT < 6 GeV/c at the same beam energy. 
This baryon- meson effect observed in p + p collisions was 
argued to be a simple reflection of the underlying dynam- 
ics of the collision. In order to have further insight in this 
matter, the spectra of resonances should be increased to 
higher momentum, which will be accomplished in STAR 
with the installation of the barrel Time-of-Flight that will 
provide extended particle identification. 

The p°, K*, A++, E*, and A* (pr) are found to be 
centrality independent. Compared to the (px) of pions, 
kaons, and anti-protons, we measure that the (pt) of 
these resonances are approximately the same as or even 
higher than the proton {px}- The resonances (px) as a 
function of centrality follow a mass ordering. 

The p°/7r-, K* / K~ , A++/p, S*/A, and A7A ratios 
are measured and we find that the p^/tt^ ratio is in- 
dependent of centrality up to the 40-80% of the inelastic 
hadronic Au+Au cross-section and it is of the same order 
of the corresponding p + p measurement. If we speculate 
there is particle re-scattering even without the presence 
of a medium for short-lived resonances. We can interpret 
these results in terms of the re-scattering/regeneration 
scenario and conclude that in both cases the regeneration 
is the dominant process. We observe that the K* /K~ ra- 
tio is the same for p + p and the most peripheral central- 
ity class in d+Au collisions. Then, it slightly decreases to 
peripheral Au+Au collisions to a suppression in central 
Au+Au collisions, showing that the re-scattering is the 
dominant process. The I]*/A and the A*/A ratios mea- 
sured in d+Au collisions are the same as those measured 
'm p + p collisions within errors, as expected since they 
are not as short-lived as the p^ , K* or 

The RdAu of the K* , and S* scale with Nbin for 
Px > 1-2 GeV/c taking into account the uncertainties in 
the normalization. We also observed that the p° RdAu for 
Px > 1-5 GeV/c is suppressed compared to the charged 
hadrons and charged pions RdAu ■ More information may 
be obtained from the RdAu measurement if it is extended 
to higher px which will be accomplished in STAR with 
the installation of the TOP. 

The measurement of these resonances in d+Au colli- 
sions will provide reference for future measurements in 
A+A collisions. 



VI. ACKNOWLEDGMENTS 

We thank the RHIC Operations Group and RCF at 
BNL, and the NERSC Center at LBNL for their sup- 
port. This work was supported in part by the HENP 
Divisions of the Office of Science of the U.S. DOE; the 
U.S. NSF; the BMBF of Germany; IN2P3, RA, RPL, 
and EMN of France; EPSRC of the United Kingdom; 
FAPESP of Brazil; the Russian Ministry of Science and 
Technology; the Ministry of Education and the NNSFC 
of China; IRP and GA of the Czech Republic, POM of 



21 



the Netherlands. DAE, DST, and CSIR of the Govern- for Scientific Research; STAA of Slovakia, and the Korea 
mcnt of India; Swiss NSF; the Polish State Committee Sci. & Eng. Foundation. 



[1] M. Cheng et ai, Phys. Rev. D 74, 054507 (2006); C. 
Bernard et al, Phys. R ev. D 75, 094505 (2007); F. 
Karsch, "hep-ph/0701210' 

[2] C. Adler et ai, Phys. Rev. C 66, 061901(R) (2002). 

[3] J. Adams et ai, Phys. Rev. Lett. 92, 092301 (2004). 

[4] J. Adams et ai, Phys. Rev. C 71, 064902 (2005). 

[5] S. Salur, J. Phys. G 32, S469 (2006). 

[6] S. Salur, PhD Thesis, Yale University 2006. 

[7] B. I. Abelev et ai, Phys. Rev. Lett. 97, 132301 (2006). 

[8] M. Bleicher et ai, Phys. Lett. B 530, 81 (2002). 

[9] M. Bleicher, Nucl. Phys. A 715, 85 (2003). 
[10] G. Torrieri et ai, Phys. Lett. B 509, 239 (2001). 
[11] J. Rafelski et ai, Phys. Rev. C 64, 054907 (2001). 
[12] J. Rafelski et al, Phys. R ev. C 65, 069902 (2002). 
[13] C. Markert et aZ., 'hep-ph/0206260 

[14] G. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 (1991). 
[15] R. Rapp, Nucl. Phys. A 725, 254 (2003). 
[16] E. V. Shuryak and G. Brown, Nucl. Phys. A 717, 322 
(2003). 

[17] J. Schaffner-Bielich, Phys. Rev. Lett. 84, 3261 (2000). 
[18] R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000). 
[19] D. Antreasyan et al, Phys. Rev. D 19, 764 (1979). 
[20] M. Anderson et al, Nucl. Instrum. Meth. A 499, 659 
(2003). 

[21] C. Adler et al, Nucl. Instrum. Meth. A 461, 337 (2001). 
[22] K. H. Ackermann et al, Nucl. Instrum. Meth. A 499, 
713 (2003). 

[23] J. Adams et al, Phy s. Rey^Lett. 91, 072304 (2003). 

[24] B. I. Abelev et ai, nucl-ex/0609021 

[25] H. Bichsel, Nucl. Instrum. Meth. A 562, 154 (2006). 

[26] Particle Data Group, J. Phys. G 33, 1 (2006). 

[27] C. Adler et ai, Phys. Lett. B 595, 143 (2004). 

[28] C. Adler et al, Phys. Rev. Lett. 89, 092301 (2002). 

[29] X. N. Wang and M. Gyulassy, Phys. Rev. D 44, 3501 

(1991); Compt. Phys. Commun. 83, 307 (1994). 
[30] C. Adler et ai, Phys. Rev. Lett. 89, 272302 (2002). 
[31] P. Braun-Munzinger (private communication). 
[32] P. F. Kolb and M. Prakash, Phys. Rev. C 67, 044902 

(2003). 

[33] C. Adler et al, Phys. Rev. Lett. 87, 112303 (2001). 
[34] M. Aguilar-Benitez et ai, Z. Phys. C 50, 405 (1991). 



[35] P. Fachini, R. S. Longacre, Z. Xu and H. Zhang, J. Phys. 
G 34, 431 (2007). 

[36] P. D. Acton et al, Z. Phys. C 56, 521 (1992). 

[37] G. D. Lafferty, Z. Phys. C 60, 659 (1993); (private com- 
munication) . 

[38] K. Ackerstaff et al, Eur. Phys. J. C 5, 411 (1998). 

[39] D. Buskuhc et ai, Z. Phys. C 69, 379 (1996). 

[40] P. Abreu et al, Phys. Lett. B 298, 236 (1993). 

[41] B. I. Abelev et ai, Phys. Rev. C 75, 064901 (2007) . 

[42] R. Arndt et ai, Phys. Rev. D 32, 1085 (1985); R. Rapp, 

Private Communication. 
[43] L. McLerran and R. Venugopalan, Phys. Rev. D 59, 

094002 (1999). 

[44] Y. V Kovchegov and L. McLerran, Phys. Rev. D 60, 
054025 (1999); Y. V Kovchegov and L. McLerran, Phys. 
Rev. D 62, 019901 (2000), Erratum. 

[45] F. Gelis and A. Peshier, Nucl. Phys. A 697, 879 (2002). 

[46] Y. V. Kovchegov and A. H. Mueller, Nucl. Phys. B 529, 
451 (1998). 

[47] A. Dumitru and L. McLerran, Nucl. Phys. A 700, 492 
(2002). 

[48] D. Kharzeev and M. Nardi, Phys. Lett. B 507, 121 
(2001). 

[49] A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D 

52, 6231 (1995). 
[50] A. Krasnitz and R. Venugopalan, Phys. Rev. Lett. 86, 

1717 (2001). 

[51] B. Alper et al, Nucl. Phys. B 87, 19 (1975). 

[52] K. Alpgard et ai, Phys. Lett. B 107, 310 (1981). 

[53] G. Gatoff and C. Y. Wong, Phys. Rev. D 46, 997 (1992). 

[54] J. Schaffner-Bielich et al. , [nucl-t h7M02054| 

[55] STAR Collaboration, in preparation. 

[56] J. Adams et al, Phys. Lett. B 637, 161 (2006). 

[57] J. Adams et ai, Phys. Lett. B 616, 8 (2005). 

[58] J. Adams et al, Phys. Rev. Lett. 92, 112301 (2004). 

[59] J. Adams et ai, Phys. Rev. C 70, 064907 (2004). 

[60] X. N. Wang, Phys. Rep. 280, 287 (1997); M. Lev and 

B. Petersson, Z. Phys. C 21, 155 (1983); T. Ochiai et al. 

Prog. Theor. Phys. 75, 288 (1986). 
[61] B. I. Abelev et al, Phys. Lett. B 655, 104 (2007). 



